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To examine the properties of the frustrated XY model at an incommensurate field we have 
examined a sequence of magnetic filling factors / which approach the irrational value of one minus the 
golden mean. At all / studied, the system undergoes a finite-temperature weak first-order transition 
involving the freezing out of fsing-like domain walls. As one approaches incommensurability, the low 
temperature phase of the system changes from the staircase states found by Halsey S to a striped 
phase consisting of a superlattice of parallel shift (Pott's-like) domain walls. 
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The frustrated XY model provides a convenient frame- 
work to study a variety of fascinating phenomena dis- 
played by numerous physical systems. One experimental 
realization of this model is in two-dimensional arrays of 
Josephson junctions and superconducting wire networks 
A perpendicular magnetic field induces a finite 
density of circulating supercurrents, or vortices, within 
the array. The interplay of two length scales - the mean 
separation of vortices and the period of the underlying 
physical array - gives rise to a wide variety of interest- 
ing physical phenomena. Many of these effects show up 
as variations in the properties of the finite-temperature 
superconducting phase transitions at different fields. Re- 
cent and ongoing experiments have been able to make 
measurements of the critical exponents in superconduct- 
ing arrays jij , opening the opportunity to do careful com- 
parison of theory and experiment. In this Letter we ex- 
amine the critical properties of the 2D XY model in the 
densely frustrated regime (/ 3> 0) and for a sequence of 
fields which approach incommensurability. 

The Hamiltonian of the frustrated XY model is 

n = -jJ2cos(e l -e 1 -A lJ ), (i) 

(v) 

where 9j is the phase on site j of a square L x L lattice 
and Aij = (2tt/4>o) f- A • dl is the integral of the vector 
potential from site i to site j with <f>o being the flux quan- 
tum. The directed sum of the around an elementary 
plaquette ^2 Aij = 27r/ where /, measured in the units of 
4>q , is the magnetic flux penetrating each plaquette due 
to the uniformly applied field. 

It is, in general, quite hard to study a system near in- 
commensurability For numerical work this is due to the 
large system sizes required; for / = p/q a system of at 
least q x q is required, and by definition q — > oo for incom- 
mensurate /. A lack of knowledge about the ground state 
of the system and the low energy excitations makes a phe- 
nomenological description of the phase transitions quite 
difficult. It has been speculated that as one approaches 
irrational /, some "soft-mode", of unknown form and 
origin, will enter the system resulting in the transition 



temperature approaching zero ||-|^] . Some works, based 
on Monte-Carlo simulations of systems near incommen- 
surability have suggested the system freezes at a finite 
temperature to a glass state MM. Our work suggests 
that this glassy behavior is due to the periodic boundary 
conditions used in these simulations. We find that, even 
near irrational /, the system has a finite temperature 
transition to an ordered state. 

To study the system near incommensurability we ex- 
amine a sequence of states, / = 3/8, 5/13, 8/21, 13/34, 
21/55, 34/89- •• which approaches the quadratic irra- 
tional value of 1 - = (3 - V5)/2 (O = 0.618 • • • is 
the golden mean). We examine ground state properties 
and low energy excitations using a numerical constrained 
optimization to minimize the energy. We correlate these 
states to those found with Monte Carlo (MC) simulations 
of systems with soft boundary conditions (soft boundary 
conditions are necessary to relax the constraint of fixed 
periodicity). We find that these systems can have low 
temperature striped phases similar to those found in com- 
mensurate-incommensurate transitions. For all / studied 
we find a finite temperature phase transition which does 
not appear to approach zero. 

The ground states of the Hamiltonian ([!]) will be 
among the solutions to the supercurrent conservation 
equations dTL/dOi = 0: 

^wa.{6j,-9 i -A i j,) = Q (2) 

j' 

where j' are the nearest neighbors to i. One set of solu- 
tions to these equations was found by Halsey || in which 
the square network is partitioned into diagonal staircases 
with a constant current flowing along each staircase. The 
resulting fluxoid patterns consist of diagonal lines of vor- 
tices. A unit cell of the staircase fluxoid pattern for 
/ = 8/21 is shown in Fig[l](a). For / = p/q, the staircase 
pattern can sit on q sub-lattices and, in addition, there 
are q states with the lines of vortices going along the op- 
posite diagonal, making a total of 2q degenerate states 
(/ = p/q 7^ 1/2). While these staircase states are not, 
as will be shown later, the ground states for all /, they 
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are a useful set of states in that the striped phases we 
find near / = fi can be defined in terms of domains of 
staircase states separated by parallel domain walls. 

Figure |l|(b) and (c) show some of the low energy do- 
main walls for a typical case, / = 8/21. Domain walls 
between the 2q degenerate staircase states can be clas- 
sified into two types. Shift walls involve a shift of the 
vortex pattern across the wall (such as in Fig ||(c) where 
the pattern on the right is shifted down by 8 lattice spac- 
ings with respect to the pattern on the left) but the lines 
of vortices are still going along the same diagonal. Her- 
ringbone walls are walls between states with the vortex 
lines going along opposite diagonals. 

To calculate energies of different vortex patterns, we 
solved equations (Q) numerically, using a quasi-Newton 
method, on lattices with up to 2.3 x 10 5 sites and con- 
straints fixing the fluxoid occupation of each plaquette. 
For f — 3/8 the lowest energy wall is a herringbone wall, 
but there is a shift wall with only slightly higher energy. 
For the higher order rationals (/ = 5/13 to 34/89) there 
is at least one shift wall with lower energy than any sin- 
gle herringbone wall. In addition, a striped phase such as 
the one shown in Fig. ||(a), consisting of a super-lattice 
of parallel shift walls is lower in energy than the plain 
staircase state for / = 8/21, 13/34, 21/55 and 34/89. 

The energy of a superlattice of parallel shift walls is 
shown in Fig ^(b) for / = 5/13 and 8/21. The presence 
of a short-range energetic repulsion makes wall crossings 
energetically unfavorable. Also, the interaction is essen- 
tially flat at large distances. There is, however, a mini- 
mum in the interaction potential at a finite separation of 
the walls. This minimum arises due the directionality of 
the wall which causes the distortion of the phase to be 
asymmetric on each side of the wall. As a result, when 
the distortion of the phases from two walls overlap, there 
can be some cancellation. If however, the walls get too 
close, the distortions of the phase field from the two walls 
start to match, causing a rapid increase in energy. 

While energy calculations show that the striped phase 
is lower in energy than the plain staircase state for 
/ = 8/21, 13/34, 21/55 and 34/89, this does not neces- 
sarily mean that there is not some other state with even 
lower energy. To test this, we undertook extensive Monte 
Carlo simulations. Boundary effects can propagate quite 
far into the system as the boundary can easily induce a 
high energy domain wall into the system. To alleviate 
this strain, the system will break this single high energy 
interface into numerous lower energy walls resulting in 
a complicated structure of walls which can extend well 
into the system. These trapped domain structures will 
have a significant effect on finite temperature states and 
transitions found in smaller samples. Free boundary con- 
ditions can also induce domain walls, as a free boundary 
can act like a mirror plane in the system. The problems of 
these boundary conditions can be alleviated by installing 
a boundary layer at the edge of the system. In the bound- 
ary layer, the coupling J in ([!]) goes from one on the in- 
terior side to zero on the exterior: J(x) = A(l — e~ x l x ) 



where A = 1/(1 - e~ wX ) and w is the number of rows 
of lattice sites in the boundary layer and A < w is ad- 
justed so that connection to the interior (where J = 1) 
is reasonably smooth (see inset on Fig. ||(a)). We found 
that w = 8 and A = 2.5 gave reasonable results. In prac- 
tice, it is only necessary to use boundary layers in one 
direction and periodic boundary conditions can be kept 
in the other. This does lead to a preferred direction in 
the striped phase (stripes parallel to the boundary layers) 
but does not seem to have a qualitative effect on the sys- 
tem behavior (it does however affect finite size effects). 
Measurements of the energy, order parameter, etc. were 
made only on the interior, where the coupling J = 1. 

For the discrete degrees of freedom we kept track of 
an orientational order parameter measuring whether 
the vortices are preferentially arranged along one diago- 
nal, and in the striped phase an order parameter p mea- 
suring the density of shift walls. The MC simulations 
used a heat bath algorithm with system sizes 32 < L < 
96. We computed about 10 7 MC steps (complete lat- 
tice updates), and data from different temperatures was 
combined and analyzed using histogram techniques ^ . 

At the lowest temperatures of the simulations, 
ksT/J — 0.03, we find the system goes into the states 
expected from the energy calculations: / = 3/8 and 5/13 
are in the plain staircase states and / = 8/21, 13/34, and 
21/55 are in a striped phase. In the following discus- 
sion, we start by examining / = 8/21, 13/34, and 21/55 
which undergo a first order phase transition at about 
ksT c — 0.13 J from the striped phase to the diagonally 
disordered phase. We then turn to the / = 5/13 case 
which has a transition from the plain staircase state to 
the striped phase at fc_eT c = 0.04J and then has a tran- 
sition at ksT c = 0.13 J to a diagonally disordered phase. 
For the largest system sizes, / = 3/8 appears to undergo 
a single transition from the plain staircase state to the 
diagonally disordered phase at about ksT c = 0.123J. 

Figure ||(a) shows the diagonal order and shift wall 
density as a function of temperature for / — 8/21. In 
the high temperature phase, domain walls of all types, 
shift and herringbone, are present and the vortex lattice 
is disordered. At the critical temperature T c , the sys- 
tem orders by freezing out herringbone walls, leaving a 
diagonally ordered striped phase. This striped phase has 
a density p = 21/29 of shift walls almost independent 
of temperature (see Fig. ||(a)). (p = 21/29 corresponds 
to an average spacing of I55-.) The non- integer spacing 
comes from a mixture of wall spacings di of 1 and 2, 
arranged in a Fibonacci sequence cut off at 21, 

d = 1, 
(di,da) = (2,l), 
(d 3 ,d 4 ,d 5 ) = (do, {d 1 ,d 2 )), 
(d 6 ,d 7 ,d 8 ,d 9 ,d 10 ) = ((rfi,rf 2 ), (d.3,d4,d 5 )), 

■■■d 21 . (3) 

Note that it is the wall spacings that repeat periodically 
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every 21 walls, but the actual vortex lattice period re- 
peats every 29 x 21 = 609 lattice constants (one pe- 
riod of the wall spacings takes 29 lattice constants for 
p = 21/29). Thus, the typical period of the ground state 
vortex lattice can be of order q 2 rather that q for / — p/q. 
The spacing observed in the Monte Carlo simulations cor- 
responds to the system sitting at the minimum of the 
energy in Fig. ||(b). We should note that this is quite 
different from the normal case studied in commensurate- 
incommensurate transitions ||Io|| where there is no mini- 
mum in the interaction potential to pin the walls. 

The transition from the striped phase to the diago- 
nally disordered phase appears to be first order. This 
is indicated by the presence of a free energy barrier at 
the transition between the ordered and disordered states 
which diverges as the system size increases fTl} ]. The 
free energy as a function of energy is obtained using 
Tl(E) = - bxP L (E) where Pl(E) is the probability dis- 
tribution for the energy generated by Monte Carlo simu- 
lation of a L x L system. Figure ||(b) shows the growth in 
this barrier as the system size increases from L = 42 to 
84 giving evidence for a first order transition. The bar- 
rier is, however, growing quite slowly so the transition is 
only weakly first order and the system sizes available are 
not large enough to apply finite size scaling to confirm 
the nature of the transition. One can do an approxi- 
mate extrapolation of the measured T C (L) to obtain T c 
= 0.1325 ± 0.0007 for / = 8/21. 

The / = 13/34 and 21/55 cases also undergo what 
appears to be a first order phase transition at around 
the same T c = 0.13 from the diagonally disordered high 
temperature phase to the striped phase. The striped 
phases for these / appears to be slightly more compli- 
cated. Like the / = 8/21 case, the stripes appear to 
be mainly composed of shift-by-eight walls with a sim- 
ilar Fibonacci sequence of spacings (with spacings of 2 
and 3). However, these higher order rationals also have 
other walls which have negative energy with respect to 
the staircase state. These additional walls also seem to 
be present in a much lower density, inter-spaced between 
the shift-by-eight walls in some quasi-periodic pattern. If 
one includes these walls, the wall density is similar to / 
= 8/21 and the vortex lattices look very similar to / = 
8/21, which is probably why they have such similar T c . 

For / = 5/13 the shift walls are the lowest energy walls, 
but the striped phase costs energy (see Fig. ||(b)). The 
striped phase can exist at finite temperature however, 
due to entropic reasons which we shall discuss below. In 
the Monte Carlo simulations we see a very similar tran- 
sition for / = 5/13 (similar T c and weak first order) to 
the one seen for / = 8/21 from the diagonally disordered 
high temperature phase to the striped phase. The wall 
density in the striped phase is fixed at about p — 13/31, 
which can be constructed from a Fibonacci sequence of 
wall spacings consisting of spacings of 2 and 3 in a man- 
ner similar to that used in Eq.([|). For L = 39, and at 
about T » 0.05 the wall spacings of 2 and 3 appear to 
switch to give a slightly lower energy state at p = 13/34. 



It is unclear however if this would be the case for larger 
systems and would require further study. These two wall 
densities, p = 13/34 and 13/31 correspond to the two 
dips within the larger minima seen in the interaction en- 
ergy shown in Figure ^(b) . At a lower temperature T p m 
0.045, the system undergoes another first order transition 
from the striped phase to the plain staircase state. 

The transition from the plain staircase state to 
the striped phase is similar to the commensurate- 
incommensurate transitions studied in the context of ad- 
sorbed films pc[ , which is a second order phase transi- 
tion. In studies of these transitions, one considers the 
free energy of a single line per unit length e s . This can 
be estimated using a simple solid-on-solid (SOS) model 
of the shift line. The energy of the line, extending from 
one side of the system to the other is 

H s {z} = a\\L + a±22\z k - Zk-i\- (4) 

k 

where mi (cj_) is the energy per unit length in the direc- 
tion parallel (perpendicular) to the wall. The heights Zk, 
take on integer values. The partition function, can be 
evaluated to give the interfacial free energy per column 
|T§ e s = Thi[e CT ii /T tanh(cr ± /(2T))]. A phase transition 
occurs when e s becomes negative. If this were the case 
here, one would see a continuous rise in the shift-wall den- 
sity. What makes this case different is the presence of the 
minimum in the wall interaction potential (Fig. ^|(b)). 

If placed in a system with other shift walls, the walls 
will experience an entropic repulsion since a wall can only 
occupy the region of space between it's neighbors. To see 
whether or not this entropic repulsion is relevant, we es- 
timate if two walls remain bound together at the minima 
of the interaction potential. This is done using a SOS 
model for two walls with a binding energy equal to the 
depth of the minima in the interaction: 

H d {A, z} = ^{(2ba + U\\5 Zk>0 ) + ba\z k - z k -i\ 

k 

+(2ba + u±S Zk>0 )Ak}. (5) 

where z k is the separation of the walls (zk > 0), and Afc 
is the number of vertical steps the two walls take in the 
same direction in the k'th column (— oo < A& < oo). 
it II and u± are the binding energies parallel and perpen- 
dicular to the wall. Summing over A& leaves the par- 
tition function in the form of a transfer matrix: Z = 
^2{z k } life I** -1 - A ground state eigenvector = 
e _MZ , where 1/p is the localization length, or typical dis- 
tance separating the lines, characterizes the bound state 
of the two lines, p = defines the unbinding transition 
at T\y. Doing this, one finds an unbinding temperature of 
ksTb/J — 0.51. Below Tf,, the entropic repulsion is in- 
sufficient to push the system out of the minimum. Above 
Th, the striped "solid" phase will melt into a phase where 
the wall density changes continuously with temperature. 
Here, however, this is preempted by the entrance of the 
diagonally disordered phase at T c — 0.13 J. 
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In order for the striped phase to be stable for / = 5/13, 
where it costs energy, there must be sufficient entropy 
from the lines wandering within the region between it's 
neighbors. The energy per line at finite temperature 
can be estimated using (^) with Zk restricted to 0, ±1 
as the minimum in the interaction energy is at a spacing 
of about 2. The free energy per line per column is then 

e b = Tln[e CT ii/ T /(l + e - 2a±/T (l + y/l + 8e 2 <W T )/2)]. 

At the point where ef, crosses zero the striped phase coex- 
ists with the plain staircase state and a first order phase 
transition occurs. Taking cry = 0.041J for the shift-by- 
eight wall at the minimum of the interaction energy and 
<t_l = 0.04 J from an average of measurements of the en- 
ergy of several kinks of differing lengths, one finds that eb 
crosses zero at Tf, = 0.04J in reasonable agreement with 
the value observed in the Monte Carlo simulations. 

In experiments 0], a finite temperature second order 
phase transition is seen at / = Q. That the transition 
occurs at finite temperature is in agreement with our re- 
sults but the continuous transition appears to disagree 
with the very weak first order phase transition found 
here. However, bond disorder which is always present 
experimentally, has been shown to wipe out any coexis- 
tence region of two phases in two dimensions making all 
transitions continuous jl3|fl. 

In conclusion, we find that all of the systems studied 
undergo a finite temperature first order transition from 
an ordered state to a diagonally disordered state. The 
transition temperature is nearly constant and shows no 
signs of approaching zero as one goes to more incommen- 
surate /. As one approaches incommensurate /, the low 
temperature state changes from the plain staircase state 
found by Halscy to the striped phase. 
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FIG. 1. Fluxoid pattern for / = sr f° r (a) a unit cell 
of the staircase state, (b) A herringbone wall, and (c) and 
shift-by-eight wall. A vortex is shown as a dark square. In 
(c), the pattern on the right is shifted down by eight from 
where it would be if it had just continued the pattern on the 
left. 
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FIG. 2. (a) A section of the striped phase for / = 8/21 
corresponding to the minimum in the energy in (b). The se- 
quence of wall spacings repeats periodically. The shading is 
only a guide for the eye. (b) Energy per unit area of a su- 
perlattice of shift-by-eight walls as a function of their average 
separation for / = A and Jr. 
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FIG. 3. (a) Diagonal order (dashed) and shift wall den- 
sity (solid) versus k B T/J for / = 8/21 and L = 42,63, and 
84. Dotted line indicates a shift wall density of 21/29. Inset: 
Couplings J in a cross-section of the system for L = 42. Data 
from the boundary layers is discarded, (b) Free energy barrier 
between ordered and disordered state for / = 8/21. 
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